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ABSTRACT 

The  propagation  of  electromagnetic  waves  in  a  region  containing 
free  electi^jns  is  discussed.  The  equations  are  linearized  by  making 
the  small  signal  approximation.  First,  it  is  shown  that  these 
linearized  equations  will  have  a  unique  solution  in  case  the  time- 
varying  current  is  parallel  to  the  velocity  if  the  tangential  com- 
ponents of  either  the  magnetic  or  electric  field  and  the  normal 
component  of  the  current  are  specified  on  the  surface  bounding  the 
region.  Second,  these  equations  are  used  to  discuss  the  propagation 
of  electromagnetic  waves  in  a  circular  cylinder  completely  filled 
with  electrons.  It  is  shown  how  to  obtain  a  convergent  expansion 
for  the  field  quantities.  This  problem  may  have  applications  to 
drift  tubes, 
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I.  IM'RODUCTION  AMD  STATEMENT  OF  FROBLEM 

The  propagation  of  sinusoidal  electromagnetic  waves  in  charge  — 
free  space,  guided  or  unbounded,  has  been  studied  in  great  detail   • 
The  purpose  of  this  rsrort  is  to  extend  the  study  to  the  guided  propagation 
of  electromagnetic  waves  in  regions  containing  charged  particles,  such  as 
an  electron  beam.  Although  the  equations  for  charge  free  space  are  linear, 
the  equations  for  electron  beams  are  nonlinear.  A  similar  study  in  the 
case  of  electron  beams  is  possible  only  upon  linearization  of  the  equations. 
The  linearization  is  accomplished  by  decomposing  each  field  quantity  into 
a  sum  of  a  tijne-independent  term  and  of  a  relatively  much  smaller  sinusoidal 
term.  After  substituting  this  into  the  original  nonlinear  equations,  it  is 
assumed  that  terms  of  higher  order  such  as  the  product  of  tvio  sinusoidal 
terms,  may  be  neglected.  The  resulting  linear  equations  govern  the  variation 
in  space  of  the  electric  and  magnetic  fields,  and,  in  addition,  the  variation 
of  those  field  quantities  which  characterize  the  electron  beam,  namely, 
electron  current  density,  velocity,  and  charge  density. 

We  shall  study  the  extensions  of  the  two  following  problems  in 

charge  free  space.  First  is  the  uniqueness  problem  for  charge-free  space, 

(2) 
It  is  well  known    that  the  fields  within  a  multiply  connected  region  are 

uniquely  determined  by  either  the  component  of  the  electric  field  or  the 
component  of  the  magnetic  field  which  is  tangent  to  the  surface  bounding 
the  region,  A  similar  uniqueness  theorem  is  proven  for  the  sinusoidal 
propagation  of  electromagnetic  waves  in  electron  beams.  The  theorem  states 
that,  under  certain  conditions,  the  sinusoidal  field  quantities  are  deter- 
mined uniquely  within  a  multiply  connected  region  by  specifying  on  the 
surfaces  bounding  the  region  the  folloi\d.ng  quantities:  the  component  of 
the  sinusoidal  electric  field  tangent  to  the  surfaces,  and  the  sinusoidal 
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electron  velocity,  or  the  sinusoidal  current  density  normal  to  the 
surfaces. 

The  second  problem  is  to  determine  the  modes  of  propagation 
of  sinusoidal  electromagnetic  waves  in  a  right  cylinder  whose  lateral 
surface  is  of  infinite  conductivity.  The  extension  of  this  problem 
to  be  studied  is  that  in  which  the  right  cylinder  is  completely  filled 
with  an  electron  beam,  whose  time-independent  velocity  is  a  function  only 
of  the  axial  distance.  It  is  also  assumed  that  the  current  density  and 
electron  velocity  possess  only  axial  coniponents.  The  determination  of 
the  modes  of  propagation  depends  upon  the  solution  of  a  system  of  two 

second-order  ordinary  differential  equations,  which  are  coupled  to  each 

(3) 
other.  These  equations  are  solved  by  the  raatrizant    method  in  the  form 

of  an  infinite  series.  ¥e  prove  the  convergence  of  this  series.  The 
necessary  boundary  conditions  are  specified  in  accordance  with  the  unique- 
ness theorem.  The  results  derived  by  this  method  may  be  applied  to  the 
design  of  microwave  electron  tubes,  for  which  a  knowledge  of  the  modes 
of  propagation  is  indispensable, 

II,  FORMULATION  OF  EQUATIONS  AND  PROOF  OF  UI^QUENESS  THEOREM 
The  basic  Maxwell  equations    in  MKS  units  are: 

^  X  F   =    -/^  Ht  (1) 

Vx  H     -     J  4-<3-^  +  e^  E^  (2) 

V.  7  4.j?^    =    0  (3) 

±L    -    -2-    Te  +  Vxh]  (U) 

7   =  /oT  (5) 
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Here,  E,  H,  J,  V,  and  p     denote  the  electric  field,  .magnetic  field, 
current  density,  electron  velocity,  and  charge  density,  respectivelyo 
6L,  /^,  cr ,   -e,  and  ra  are  the  dielectric  constant,  permeability, 
conductivity  of  the  medium,  charge  and  mass  of  the  electron,  respectively, 
which  are  assumed  to  be  constants  and  positive  real  numbers.  The  sub- 
script t  denotes  partial  differentiation  with  respect  to  time,  and  — 

dt 

denotes  the  total  derivative.  Thus 

^  =  -^  +  V  .  (6) 

dt     5t 

Coital  letters  ^d-ll  denote  vectors  in  this  section,  unless  it  is 

specified  otherwise. 

In  accordance  with  the  linearization  process,  each  field 

quantity  is  written  as  the  sum  of  a  time-independent  terra  and  a  relatively 

much  analler  sinusoidal  terra.     Thus,   for  exajiple, 

1 
7    =     Jg  +  J  e:q3   (i^t)    ;   i     =     (-  1)^  (7) 

|Jl«|Jo| 

Here   <^  is  the  angular  frequency  which  is  a  positive  real  number.     The 
exponential  function  is  introduced  to  simplify  the  resultant  linear 
equations.     The  physical  solution  is  obtained  by  taking  the  real  part  of 
the  right-hand  side  of  equation   (7),     Substituting  formulas  similar  to   (7) 
into  the  original  equations  and  neglecting  higher  order  terras,  we  obtain 
the  following  equations  after  separating  the  time-independent  and  the  tine- 
dependent  terras, 

A,     Time-Dependent  Equations 

rx  E     =     -  i  O)/^  H  (la) 

Vx  E     m     J  +   <rE  +     iiOfco   E  (2a) 


-  u  - 


f'  •  J  +  i  OJ^  =  0  (3a) 


i  O)?  4  V-  '^  7  i  Y  •    V  Y 


o  •  T  •    ■   '0 


J 

= 

Po  V  4     /°V, 

Time 

-Independent  Equations 

Vx  Eg 

s 

0 

VxH^ 

■ 

Jo*-     ^^0 

^-   Jo 

B 

0 

^0 

•     ^  V 

0 

■ 

i  K* 

Jo 

= 

'  0    "o 

J  [e  +  V  X  Ho  4  Vq  X  h]  (Ua) 

(5a) 


(lb) 

(2b) 

(3b) 

X  H^^  (Ub) 

(5b) 

It  is  assumed  that  the  time-independent  field  quantities  have 
already  been  obtained  from  equations  (lb)  -  (5b)  and  hence  are  to  be 
considered  as  known  quantities.  The  problem  is  to  solve  equations  (la)  - 
($a)  for  the  sinusoidal  field  quantities  subject  to  boundary  conditions 
which  are  chosen  in  accordance  with  a  uniqueness  theorem.  In  view  of 
the  complexity  of  the  equations,  it  has  not  been  possible  to  obtain  a 
uniqueness  theorem  for  all  cases.  However  for  the  fecial  case  in  which 
J  and  V  are  parallel,  the  following  uniqueness  theorem  is  proved. 

Theorem  I 

If  J  and  V  are  parallel  within  a  multiply  connected  region  G, 
the  sinusoidal  field  quantities  which  are  bounded  and  possess  continuous 
first  derivatives  are  determined  uniquely  within  G  by  specification  of 
the  following  quantities  on  the  boundary  surfaces  3:  E  x  n  or  H  x  n, 
and  V  or  J  •  n  ,  n  is  the  unit  normal  to  3  pointing  away  from  the  region  G, 
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Proof 

From  equations  (la),  (2a)  and  the  use  of  Gauss'  integral 
theorem,  it  follows  that 

(H*  •  V  X  E  -  E  •  V  X  H*)  d  T  =  )V7-  (E  X  H*')  dT  =   E  x  H*  •  d  s 

(  U,  "  '  "*  ^"^  -  i^  ro  E  «  E*  dY  +  6'e  •  E*  dT  +  U*  •  E  dT  .  (6) 


B>  itU 


Here  d  T  is  an  element  of  volume  within  G  and  the  asterisk  denotes  the 
complex  conjugate.  Using  the  vector  equation  identity, 

y  (Vq  .V)  =  Vq  'VV  I  V  'VVq  4.  Vq  X  57x  V  I  V  X  Vx  7q  ,   (7) 

equation  (Ua)  may  be  rewritten  as  follows: 

i  u;  V  I  V(Vo  •  V)  m     -S-|^  E  +  V  X  Hq  +  Vq  X  HJ  -  Vq  X  V  X  V  -  V  X  V^  X  Vg  .   (8) 

Also  since  J  and  V  are  parallel,  J,  V,   J^,   and  V^  are  parallel  to  each 
other.     Hence  taking  the  scalar  product  of  equation  (8)   with  J  , 

J*  o  E     =     J^f  iwV   •   J*  +  J*"  «  V  (Vq   •  V)  J  .  (9) 

Using  equations  (3a)   and   (5a), 

iWV  •   J*     =     itU^o  V   »  ?*  +  V«    (Vq   f  V)    J*  -  J*  •  7(V^   .  V)o  (10) 

Substituting  (10)  into  equation  (6),  one  obtains 

iu»; te H  •  H*  -  60  E  •  E*  4-  -|-  /^o  ^  °  ^  J  '^^  ^r^ '  ^*  '^^ 

+  I  [e  X  H*  +  (Vq  •  V)  /]  •  d  S  =  0  .  (11) 

Equation  (ll)  may  now  be  used  to  prove  the  uniqueness  theorem.  Let  us 
assume  that  there  are  two  field  solutions  satisfying  the  same  boundary 
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conditions  mentioned  in  the  theorem.  The  difference  of  these  two  solutions, 
which  is  also  a  solution  because  of  the  linearity  of  the  equations,  will 
have  corresponding  boundary  values  equal  to  0,  Hence  for  the  difference 
solution,  the  surface  integral  in  equation  (11)  will  vanish.  Equating 
real  and  imaginary  parts  in  equation  (11) ,  it  follows  that 


J  <r  E  •  E'""  dr   =  0,  (12) 


which  is  possible  only  if  E  s  0,  since  the  integrand  is  positive  definite. 
The  vanishing  of  the  other  field  components  follows  directly  from  the 
linearized  equations.  That  the  solution  is  unique  is  an  immediate  conse- 
quence of  the  vanishing  of  a  solution  with  zero  boundary  conditions. 
In  the  case  of  ^T  =  0,  it  can  be  deduced  only  that 

Jf^o  H  •  IT"  -  e^  E  .  E""  4-  21  /^^  V  .  V"']  dT  =  0  .   (13) 

This  equation  can  be  satisfied  with  non-vanishing  fields  only  at  certain 
discrete  frequencies,  which  are  the  resonant  frequencies  of  the  region  G 
in  the  presence  of  an  electron  beam.  In  this  case,  the  uniqueness  theorem 
will  hold  at  all  frequencies  except  those  which  are  resonant. 

The  uniqueness  theorem  has  been  proven  only  for  field  solutions 
in  which  J  is  parallel  to  V,  An  important  question  is  whether  there 
exist  non-trivial  cases  in  which  the  assujnption  prevails.  This  case 
occurs  when  H^  is  infinite.  Thus,  if  it  is  assumed  that  Vq  is  bounded 
and  possesses  bounded  derivatives,  it  follows  from  equations  (Ua)  and  (Ub) 
that  both  V  and  Vq  are  parallel  to  Hq.  Hence  J  is  parallel  to  V,  This 
case  is  discussed  in  greater  detail  subsequently. 
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III.  MODES  OF  PROPAGATION  IN  A  RIGHT  CILINDSR  COMPLETELY 

FILLED  V.TTH  .4N  ELECTRON  BEm. 

The  problem  to  be  considered  is  to  obtain  solutions  of  the 
linearized  equations  within  a  right  cylinder  of  arbitrary  cross  section 
between  the  planes  z  a  z  ,  z-^  (Fig,  1). 

I 
I 

I 

i 


*Ii.  — '^  o«'S 


The  lateral  surface  is  assumed  to  be  of  infinite  conductivity 
and  Hq  is  assumed  to  be  infinite  and  along  the  z  axis.     On  the  planes 
z  =  Zq,   z-^,  E  X  £  and  J  are  given;   z  being  the  unit  vector  along  the  z 
axis.     Another  assumption  is  that  V^  is  an  analyiiic  function  of  z  only 
and  does  not  vanish  anywhere  in  the  interval  z ^  z  /.  Z]_,     As  shown 
previously,   J,  V,  Vg,   and  J^  are  parallel  to  the  Z  axis.     Thus  the 
boundary  conditions  in  the  region  G  are: 

at  z     ■     Zq      :     E  X  z_  and  J     =     Eg   (0)    and  J  (0),  respectively;      (ih) 

z     ■     Z-,      :     E  X  z_  and  J     «     Eg   (l)   and  J  (l),  respectively;      (15) 

on  the  lateral  surface:  Exn»0  J«n=0,  (l6) 

Here  n  denotes  the  unit  vector  normal  to  the  boundary  curve  of  the  cross 
section  and  the  subscript  s  denotes  confjonents  of  the  vector  in  the  trans- 
verse plane. 
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It  follovs  from  the  uniqueness  theorem  that  there  exists  at  most 
in  the  region  G  one  solution  only  «hich  satisfies  these  boundary  conditions. 
Therefore,  any  solution  which  is  constructed  with  these  boundary  conditions 
must  be  the  only  one,  A  solution  will  now  be  given. 

The  following  additional  notation  will  be  used.  Let  J,  V,  V^,  J  , 
E,  and  H  denote  the  z  conqjonents  of  the  respective  vectors.  Also  the  sub- 
script z  will  denote  differentiation  with  respect  to  z.  Then  it  follows 
from  equations  (3a)  -  (5a)  and  from  the  fact  that  Jq  is  constant, 

P--    -      i  J,  ,  (17) 

o  -• 

V^3  j^^  ^  [2  i«  V^3  ^  3  ^2    v^  J  j^  +  [2  i  CO  V,  V^^  -  ^^  V,] J 

s   _o    i  CO  E  .  (19) 

m 

From  equations  (la)   -   (2a),  one  obtains 

\z+^s'    E  +  k^E     =       J^       [k2  J4  J  (20) 

^s  •  ^s    =    -    ite      ^z  -  \  (21) 

k^  ^s  V  ^szz     =     (V3  E)^  f  iu>Ao   (£  X  ^3  H)  (22) 

^s     '     '      r~Z    [(^^  Vz-^^V^e]  (23) 

•  (0 

H       +  V^^    H  +  k^  H     =     0   .  (2li) 

zz         s 


Here,  e  =    e     _  i  -2-.  . 

O  i4> 


CO 

k2  =  a>2  eu   ,  (25) 


and   V       =     gradient  operator  in  transverse  plane. 
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With  these  equations,  all  the  field  quantities  are  expressea 
in  terms  of  J,  E,  H,  and  their  derivatives.  It  may  also  be  shovn  by 
reversing  the  algebraic  steps  in  the  derivation  that  equations  (17)  - 
(25)  are  equivalent  to  the  linearized  equations  if  the  field  quantities 
possess  continuous  second  derivatives, 

(U) 

It  is  convenient  to  split  Eg  and  H  into  two  parts    depending 
upon  E  and  H,  respectively.  These  are  known  as  E  modes  and  H  modes, 
respectively.  Let 

E3  =  E^^^  ^e}^^    ,  (26) 

and  similarly  for  Hg,  These  quantities  must  satisfy  the  following 
equations : 

^3  •  ^s^^  =  -   T^  ^z'\  (27) 

^'  ^i^^  ^  ^sS  -  (^s  ^)z  (28) 

"i"^  -  -   iTT:  f  (£  -  ^3^'\  -  i  -  V3  e]  (29) 

i?   E^®  t  E^<^'  .  i  ul^^   (z  X  Vs  H)  (30) 

(^)  ,  _   T^  [  -  eJ2']  (31) 


^S 


Writing  each  field  quantity  as  a  product  of  two  factors  depending 
upon  coordinates  in  the  transverse  plane  or  z  only:  e,g, 

E  =  E(s)  E(z)  (32) 

« 
the  equations  separate  into  two  sets  governing  the  dependence  upon  s  and  z 

respectively.  These  are: 
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A.  E  Mode 


1,  Dependence  Upon  s 

2      2 
7   E  ♦  T,   E  =  0 
s      1 

J  =  E 

(33) 
Et  =  V3  E 

Ht  =  z  X  V3  E 

E  s  0  on  boundary  curve  of  cross  section.  Here  the  variable  6  is 
understood, 

2,  Dependence  Upon  z 

eJ 
«     _jo  i  ^  E  (3U) 

m 

\z  +  (k^  "  -!  )  E    =    ^  (k'  J  ♦  J,,)  (35) 

_T^2     ,a)    c     -     ^^        J,-E^  (36) 


To 


Here  the  variable  z  is  understood, 
B,     H  Mode 

1,     Dependence  Upon  s 

^g     H  4  T^  H     s     0 


E  (2)     .     z  X  V     H  (38) 

s  ^         s 


H 


and  £  X  Vg  H     •     0  on  boundary  curve  of  cross  section. 
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2,  Dependence  Upon  z 

H  ^  +  (k^  _  T  2)  H  =  0  {39) 

H  =  exp  (£  I  h^  z)  J  h^^  =  k^  -  T2^        (UO) 

k  -  h2 

The  boundary  conditions  in  equations  (33)  and  (38)  ensure  that 
the  boundary  conditions  in  equation  (16)  are  satisfied  for  the  component  of 
electric  field  tangent  to  the  lateral  surface.  Thus  the  only  remaining 
boundary  conditions  to  be  satisfied  are  those  on  z  •  ^o*  ^n 

Now  it  is  well  known  that  the  E(s)  and  H(s),  whose  variations 
are  given  by  equations  (33)  and  (38)  \-ilth  their  respective  boundary  conditions, 

<*>  0  0 

exist  onlj*-  for  a  discrete  set  of  values  of  T-^    and  T2  ,  namely  T-j^^^  and  T2JJ, 

(5) 
which  are  real  positive  numbers   .  Furthermore,  the  eigenfunctions  of  E(s) 

and  H(s)  corresponding  to  different  eigenvalues  are  orthogonal  to  each  other 
and  there  exists  an  orthogonal  subset  of  the  eigenfunctions  in  terms  of 
which  any  function  of  s  possesses  a  unique  convergent  expansion   •  Similarly 
the  corresponding  E^  (s)  and  E^^  -^(s),  and  H^^  \s)   and  H^^^^  (s)  form  a  complete 
orthogonal  set.  These  statements  can  easily  be  shown  to  be  true  in  the  case 
where  the  cross  section  is  a  rectangle  or  a  circle  by  using  trigonometric 
or  Bessel  eigenfunctions,  respectively.  It  vail  be  assumed  that  for  other 
cross  sections  which  are  simply  connected,  it  is  still  possible  to  expand 
Eg(l),  Eg(2),  and  J(l),  and  J(2)  in  terms  of  transverse  and  axial  eigen- 
functions respectively  and  treat  each  term  in  the  expansion  separately. 
Expanding  the  fields  in  the  region  G  in  terms  of  these  eigenfunctions,  the 
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coefficients  will  be  functions  of  z  whose  variation  is  given  by  the  Diode 
equations  with  T^^  and  1^   replaced  by  T-j_^  and  '^^n*     '^^^     ^^^* 


E  =  2  E„(s)  E„(z) 


f7 


(U3) 


Es  -  2  E3^)(s)E3(%)+2Ei^(s)E^2)(,) 
Hs  =  2  H3„^^  (s)  ^P   (z)  4ZHi2)(3)  4^)i,) 

E,(l)  =  2  E^^Cs)  a^(l)  +  2  ^^^\s)   a^(l) 
n  m 

J(l)  =  £  E^(s)  b^(l) 

and  similarly  with  Eg(2)  and  J(2),  The  a^^,  aj„,  and  b^  are  constantso 

The  coefficients  for  the  H  mode  are  easily  obtained  from  the 
H  mode  expressions  in  equations  09)  -   (U2) .  It  is  seen  that  there  are 
two  values  corresponding  to  h2  «=  -f  h2m  to  be  determined.  These  are  linear 
combinations  of  a^(l)  and  a  (2).  They  could  also  be  determined  if  both 
E_  and  H  are  specified  in  z  ■  z  and  nothing  is  specified  at  z  =  z. 

The  coefficients  for  the  E  modes  are  to  be  obtained  from  the 
expressions  for  the  E  mode  in  equations  (3U)  -  (37)  with  the  subscript  n 
added  for  the  respective  mode.  Deleting  the  subscript  n,  these  constitute 
a  set  of  ordinary  differential  equations,  iji  which  the  values  of  Eg(z)  and 
J(z)  are  specified  at  z  =  Zq,  Zn .  This  is  equivalent  to  solving  the  set 
of  ordinary  homogeneous  differential  equation  in  J  and  E  subject  to  the 

boundary  condition  that  the  values  of  J  and  J^  *.  1  o>  €.  E^   are  given  at 

(7) 
z  =  z  ,  z, .  This  constitutes  a  linear  differential  system^''  with  a  two- 
0    i 

point  boundary  condition.  The  necessary  and  sufficient  condition  that 
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there  exist  a  unique  solution  is  that  the  homogeneous  differential  equations 
Tvith  the  boundary  conditions  J  and  J^lw^S  =Oatzs  Zq,  Z]_  possess 
only  the  trivial  solution  J  s  E  «  0,  This  has  to  be  the  case;  otherwise, 
there  would  exist  a  solution  E(s)  J(z),  E(s)  E(z)  in  region  G  (Fig,  1) 
for  which  Eg  and  J  are  zero  on  the  boundary,  wliich  is  contrary  to  the 
uniqueness  theorem.  Thus  the  solution  of  the  inhomogeneous  boundary  value 
problem  exists  and  is  unique.  The  folloxving  theorem  may  now  be  proved. 

Theorem  II 

The  fields  in  G  exist  and  are  uniquely  determined  by  the 
boundary  condition  that  the  value  of  J,  J  or  V,  Eg  and  H  be  given  at 
z  =  Zq.  Nothing  need  be  specified  at  z  =  z^ . 
Proof 

It  is  obvious  from  the  expressions  for  the  E  mode  in  equations 
(3l()  »  (37)  that  J,  J  ,  E  and  H  at  z  »  z  determine  uniquely  the  values 
of  J,  Jj,,  E,  and  E_  at  z  •  z  and  conversely.  Then  according  to  the 
existence  theorem  of  ordinary  differential  equations  with  analytic  coefficients, 
there  exists  a  unique  solution  which  determines  the  value  of  J  and  Eg  at  z  ■  z, , 
Hence  the  field  solution  is  unique  according  to  the  existence  theorem  1.  The 
following  corollaiy  is  easily  proved. 
Corollary  1 

The  fields  in  region  G  e>d.st  and  are  determined  uniquely  by  the 
following  boundary  conditions: 

a)  either  H  or  E„  and  either  J  or  V  given  at  z  s  z^,  z-. 

'  S      S  O'   X 

b)  J,  J  or  V,  E,  and  E  given  at  z  s  z 

c)  J,  J  or  V,  E  and  Hg  given  at  z  -  z 
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The  remainder  of  the  thesis  vill  now  be  devoted  to  a  study  of  the  solution 
of  the  ordinary  differential  equations  for  J  and  E« 

IV.     SOLUTION  OF  THE  DIFFERENTIAL  EQUATIONS  FOR  J  AND  E    ' 

It  is  convenient  to  -write  the  ordinary  differential  equations 
for  J(z)   and  E(z)   in  matrix  form. 


U. 


Here  U  is  the  column  matrix 


U 


}    £  =  i     e  E 


(UU) 
(U5) 


,2 

'0 


T    +  i  2T 

'n  oz 


0 


-  i  2r+  e 


T, 


oz 


0         T 


.k2-T^-i2y    ,       i2r-3^ 


2     2 
^    ^o 


2        2       2   2 


(1;6) 


K'- 


V. 


2     _ 


O)' 


0 


(U7) 


The  boundary  condition  is    U     ■    U(Zq)   at  z     =     z^.  (U8) 

A  simple  solution  is  obtained  only  in  the  case  Vq     b     constant  and 

/ON 

iy^l  ^^  1   .  In  this  case  the  solution  may  be  expressed  as  the  sum 
of  four  exponential  functions  with  constant  coefficients  as  follows : 


E(z)  =   y-   E  er-piX   z),  m  =  1,  2.  3.  ^ 

I —    m     m 
m 

(-V^fk2-T|) 


(U9) 


(^0) 
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2  ^2 

^3,  U  =  *  i  h.  i-  —2 1 ..   -i—  y2  (51) 

~         (£  ^1  +  l^o)      ^  \ 

h-L  s  (k^  -T^2)*  .  (52) 

I  2 1       ^oz 

We  consider  the  problen  in  which  |  y  |  and  — —  are  both  small.  The 

o 
previous  case  suggests  finding  a  solution  expressed  as  the  sum  of  four 

exponential  functions,  whose  coefficients  are  power  series  in  y  and  V  • 

A  convenient  method  to  accomplish  this  is  by  the  matrizant 

scheme  of  Peano  and  Baker   '    and  its  extension  by  Keller  and  Keller^  ' • 

According  to  this  method,  the  solution  of  the  differential  equation  {hh)   is 

z 
u    =    ^   {  a]  U(z^)  (53) 

where 

^      i    }  /?  r?  /^2 

-0- U     =   matrizant   of  A  =   1  +   rA(x^)dx^  +yA(i2)di2/       A(x^)(ix^  + (5^*) 

°  0  0  0 

Inasmuch  as  VI (z)  ip  an  analytic  function  in  the  region 
Zq  ^  z  ^  z-L,i2|A  I  exists  and  converges  uniformly  in  the  region.  The 
following  properties  of  the  matrizant  will  be  used  here: 

a)  n      -{aI  exists  if  i^{A|  exists^  .  (55) 

b)  For  any  non^singular  matrix  P, 


A  {aJ=  P  a    |p-'  AP  .  T-\   \     P-^  (zj^^°^ 


(56) 


^0 


c)     For  any  diagonal  matrix  A  «  diag,   (J^i ,  <^ 2  >  ^^3 f '^ U^ ' 
0  00 


(10) 


(57) 
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^  (9) 

Zq     *  '  Zo    C       J    Zq     I     Zq       ^     '  Zo    «>       ■'J 


(58) 


Using  equation  (^6),  the  solution  is 


2      f   -1 
U     =     P  XI    f  P 
z„   ^ 


A  P  -  P~     P, 


]     P-^(^o)  ^(^o^    • 


(59) 


Usually  P  is  chosen  to  diagonalize  A.  Hovever,  in  order  to  separate  out 
the  effects  due  to  Vq^j  let 


A  =  A 


1  +  ^: 


(60) 


h 


0 

1 

0 

0 

< 

-12  T 

0 

/i 

0 

0 

0 

0 

1 

-^^-rl 

12  Y^ 

2     2     2     2 

0 

(61) 


0 

0 

0 

0 

12  T 
o 

31ogT^ 

0 

0 

0 

0 

0 

0 

i2T„ 

-31ogT„ 

0 

0 

(62) 


Then  by  equation  (60) 


U  =  P   XL  '{p-^P  +  p-^A^  p  _  p-lp 


^rhz^MT^)       (63) 
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where  P  is  chosen  to  diagonalize  A]_.     Such  a  non-singular  matrix  P  exists 
if  the  eigenvalues  of  A|   are  distinct. 

The  eigenvalues  of  A.-^,    A-  ,   satisfy  the  fourth-degree  equation, 


,2 

t 

2 


2  U^  +  k2  .T{)    (A+  ±Y^f 


The  roots  of  this  equation  when  \  J  \     is  small  are  obtained 

%  (0) 
as  follows:  for  y  «  0,  the  roots  K^         are: 


^°^     «  -iVo,  -iVo,  i  V  -ih-i 


^1,2,3, 


hi  -  (k^  -Ti^)^ 

No\j   \j^     is  not  zero  atA.=/\o  »  ^  ]        •     Hence,  according  to  a 

well  known  theorem  in  the  theory  of  functions  of  a  complex  variable,  \ 

is  an  analytic  function  in  the  neighborhood  of  y^  ■  9  and  the  corresponding 


power  series  is 


;\  «  X  f,^  +  /  1 +  0  (y^)  .  (65) 


^'^  ^ '  ^) 


"dK" 


_       ^  T  ^ 
Hence,  X     =  *•  i  h,  4  — ^2-_L-_-    .-i-   y%  0  (y^)  .  (66) 

3,u         -  (fhi4ro)     2h3_ 

The  remaining  two  roots  are  obtained  by  extracting  the  square  root  of 
equation  (U3)  and  applying  the  same  reasoning, 

(X^  +  ic^  -  T^)^/2(^  ^  .^^) 
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and  \ 


1,2 


=  -  1 


y  1  Y     2 

'^  o  I.  '  o  ^^ 


y  +  O(y^)  , 


(67) 


/  ^  2   ,2   „  2s  2 

(-  r/   +  k  -  T  ) 

0         i 


Thus  it  is  seen  that  there  is  a  region  in  the  neighborhood  of  y  -  0, 
excluding  y  =  0,  in  which  the  roots  are  distinct.  Also  let 


(^l,2*i^o) 


2  _ 


=1,2  '  o 


vr  2   2 


^1,2 


2    2 


r-o  4k^ 


+  V^^^' 


-  T, 


(68) 


The  matrrx  P  is  obtained  as  follows:  The  elements  of  P  or  P   satisfy 

mn      '' 

the  equations. 


In 


2  2 


2n 


Jn. 


yV''^   (X  +ir  )^  X  (x  +ir  )^  ' 


o  n 


n  n 


n  =  1,2,3,'^. 


(69) 


Hence, 


P  = 


1 

1 

2     2 
y    ^0 

\ 

^2 

2   2 
^•^0^3 

^1 

g2 

(X3.iT^) 

^^1 

Xggg 

H^3^^V 

2  2 


(70) 


To  obtain  P  ,  let  P  =  P^^  |  y^  Y^^^   P. 


(71) 


1 


\Zj 


\S2 


0 
0 


(X3.iY^)' 


0 
0 


(72) 
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0 

0 

1 

1 

2= 

0 

0 

s 

X 

0 

0 

0 

0 

0 

0 

0 

0 

(73) 


It  is   convenient  to  write  these  matrices  in  terms  of  2  x  2  submatrices 
as  follows: 


"ll         ° 


n^.       "23 


^2  = 


0       n 


21 


0       0 


(7ii) 


,-1 


-1 

^11 


-1  -1  -1 

■"13  "12  "11  '  "l3 


(75) 


and 


'11 


1        1 


\ 


12"  "11 


^1     ° 
0     go 


;  ^^13= 


1     1 


^^3    \ 


(X^^irj 


\^^V' 


(76) 


-1        1 


'11      \^-X^ 


X^     -1 


'\ 


-1 

^12 


-1 


0         g, 


-1 


(77) 


-1        1 


'13      Xi^-X3 


(X3.iT^)' 
0 


0 


(V'^o'' 


-X. 


-1 
1 


(73) 


'21 


1        1 


^3      \ 


-1        1 


'       '21       X|  -X- 
h    3 


-h 


1 
1 


(79) 
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Hence, 

P  =  Pi[i4/r,2  P^-^  pj 
p-^  =  Li  *  7\y,'^,T'  PI'  =  PI'  ^  ^  ^^'^'  (80) 

where      Hv^)   =  ^  (-D^  (j't^)"  ^h\^'^l' 

The  expansion  is  valid  and  uniformly  convergent  as  long  as 

I  ^  y^  >"o^  M  I  <  1  (81) 

vhere  M  is  the  least  upper  bormd  of  the  elements  of  P-j_  P  ♦  Note  that 
the  elements  of  ?■[  and  Pp  are  bounded  for  all  values  of  y  except  that 
of  ff    ~   ,  The  elements  of  "^y^      are     0(-— -),  Hence, 

M  =  O(-l-)  (82) 

and       I  U  y^  Y^^     wj  =  0(y) 

Thus  by  taking  y  anall  enough,  the  inequality  (81)   will  be  satisfied. 

The  matrix  R  of  which  the  raatrizant  is  sought  in  equation  (63) 
may  novj  be  expanded  in  powers  of  y^  as  follows: 

A=     P"-^  A.    P     =     diag.    (A-,,  X       ^^,A>) 

i  J.       d       J       u  (33) 

^=     P"'  ^2x  P  =     Pl""^  A2x  Pi  44^  (y^  V„^)    , 

where  H^  (y^  V^^)  «    E   {-if  {j^  vX    (Pf^  P^^'^L  Pi"'  ^2  ^-^  hx  h 

-^l"^^2x^l      ''        (^^^ 
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where     9 


^  „       p       9  n        _i         n-1  r    -1  »i  Y^^        -1 

=  2  (-1)"  (/  r^)   (p^    pp     \_\    p^  Pi    Pi,  -  3 12-.  p^    P^ 

Thus, 

R     =A+  R^  4  R^     =^4  P^"^    A2^  Pi  -  P^-^  Pi^  4t(y^  V^^)  -  9  (y^,  y^^) .   (86) 

(10)    • 
It  has  been  shown  that  it  is  possible  to  express  the  solution 

of  the  differential  equation  as  a  sum  of  e:^onential  functions  bj  letting 

R     =     D  4  R3  (87) 

D  is  a  diagonal  matrix     ~     diag,   (d  ) 
Then  using  equations (57)   and   (58), 

(88) 


^[r|    =    .a(D|i2[M| 


[dI    =     diag./exp      /        dj„  d  X  I,     m  =  1,2,3, U.     (89) 


Now  let  the  elements  of  these  matrices  be  chosen  as  follows: 


d       =  Ia  i  P-,'"^   A-    P,  -  P  "-^  P,     I  (90) 

run        [-^^  ^    1        2x    1        1        Ix  J  ^ 

(Ro)         =     (1  -  (J     )    fp,"^   A^     P,   -  P  "^  P.  1         ♦•'t'„„  -  9„„       (91) 
■>  xm  Din^    (_l         2x1         1         IxJinn^mn        mn"' 

^0 
Also  let    -t2|  m]    =     1  +    it 

i  ^     -     f        V  ^  ^  +    ^    J         \k  (^)   d  X2  J     ^  Mk„  (x^)   dxi  (93) 

Zo  °  ^0  
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4^    is  the  sura  of  tvo  terms,  of  liiich  one  is  a  combination  of  the 
J-  ran  ' 

2 
derivatives  of  V  and  the  other  is  a  power  series  in  y  ,  The  solution 

0 


of  the  differential  equation  is  now  given  by 

:9U) 


U  o  P   |d  ■  [Pi^(Zo)+  i   Pi'^(z^)*  ^(^o^^  ^  *^^^o^  "^^o^*   ^ 


The  elements  of  U  are 


V  E  fc  (fl'(-o'*4!fl'(-o'*  ♦(-o)*i*<^o»jk  \^-o^'„f^/^i^ 


(95) 


Thus  the  following  theorem  has  been  proved. 

Theorem  III 

The  solution  of  the  differential  equation  (Uh)   may  be  e:q)ressed  as 
the  sum  of  four  e3q)onential  functions  with  coefficients  given  in  the  form 

of  uniformly  absolutely  convergent  series  in  the  region  z  <  z  ^  z,  in 

2 
powers  of  y  and  V  • 

V.  SUMMARY  OF  RESULTS 

The  uniqueness  theorem  for  the  electromagnetic  field  quantities 
in  charge  free  space  has  been  extended  under  certain  conditions  to  the 
field  quantities  in  an  electron  beam.  The  uniqueness  theorem  states  that 
if  J  and  V  are  parallel  within  a  multiply  connected  region  G,  the  sinusoidal 
field  quantities  which  are  bounded  and  possess  continuous  first  derivatives 
are  determined  uniquely  within  G  by  specification  on  the  boundary  surfaces 
of  the  component  of  E  or  the  component  of  H  tangent  to  the  boundary  surfaces 
and  the  velocity  or  the  conponent  of  J  normal  to  the  boundary  surfaces. 
This  uniqueness  theorem  is  utilized  to  obtain  the  boundary  conditions  which 
determine  a  unique  solution  for  the  field  quantities  in  a  right  cylinder 
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completed  filled  vdth  an  electron  beam.  The  theorem  for  the  required 
boundary  conditions  states  that  if"  V  is  an  analytic  function  of  the 
axial  dimension  only  and  Hq  =  ao  ,  the  fields  within  the  right  cylinder 
exist  and  are  uniquely  determined  by  the  boundaiy  condition  that  the 
value  of  J,  Jg  or  7,  and  the  components  of  E  and  H  in  the  transverse 
plane  be  given  on  a  transverse  plane,  J  and  the  axial  component  of  E 
are  determined  by  a  system  of  tvo  second-order  ordinary  differential 
equations,  which  are  coupled  to  each  other.  This  system  of  differential 
equations  is  solved  by  the  matrizant  method,  and  the  solution  is  expressed 
as  the  sum  of  four  exponential  functions  with  coefficients  given  in  the 
form  of  uniformly  absolutely  convergent  series  in  powers  of  y^  and  Vq2, 

The  author  wishes  to  thank  Professor  B,  Friedman  for  his 
constant  encouragement  and  aid  in  carrying  out  this  work. 
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